arXiv:1505.07617vl [physics.optics] 28 May 2015 


Isolation of coherent and incoherent nonlinear spectroscopic 

signals by phase modulation 

Khadga Jung Karki/’[^ Loni Kringle,^ Andrew H. Marcus,^ and Tonu Pullerits^ 

^Chemical Physics, Lund University, Getingevdgen 60, 222 4 I, Lund, Sweden 
‘^Department of Chemistry and Biochemistry, 

Univerisity of Oregon, Eugene, OR 97403-1253, US 

We investigate the effect of phase modulation of laser beams on the coherent 
and incoherent non-linear signals arising from the interaction of femtosecond pulses 
with matter. We observe that two collinear beams, whose phases are modulated by 
frequencies fi and (/> 2 , produce two second harmonic signals from non-linear crys¬ 
tals whose intensities at the detector are modulated at the frequencies (/>2 — 4>i and 
2{(j)2 — 4>i)- We also observe that an incoherent action signal, such as fluorescence and 
photocurrent, which arises from the absorption of two photons, is modulated at the 
same frequencies as in the case of second harmonic generation. We present a theoret¬ 
ical analysis to explain our observations. These results are important to understand 
how phase modulation techniques can be used to isolate different field-matter inter¬ 
action pathways in a non-linear process. Because the method uses modulation of 
the signal intensity rather than wave-vector matching to isolate different signals, it 
could be useful to perform multi-photon absorption studies on single molecules or 
nanoparticles. 
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I. INTRODUCTION 

Pulsed lasers, which are routinely used in spectroscopy, can have very high peak power. 
For example, the peak power of a 10 fs pulse with an energy of 1 nJ is on the order of 100 
kW. When such pulses are focused to a spot size of 10 pm the peak intensity can reach ~ 100 
GW/cm^. Interaction of pulses having such high intensity electric helds with a sample can 
generate a plethora of non-linear signals, which can be used to study various spectroscopic 
properties of the sample. BIZI It is usually cumbersome, and also unnecessary, to investigate 
all of the non-linear signals in an experiment. Hence, different methods, such as phase 
matching, phase cycling and phase modulation, are often used to isolate the relevant signal 
contributions. 

For an electric held described by a plane wave E(t) = A(t)[exp{i{kx — ut)) +exp{—i{kx — 
cut))], where A{t) is the temporal envelope of the held, the wave-vector k and the frequency 
uj quantify the momentum and energy of the photons, respectively. Conservation of energy 
and momentum in an u-wave mixing process |6] leads to the generation of phase-matched 
non-linear signals with distinct wave-vectors kj = ±ki ± ^2 ± ... ± fc„_i and frequencies 
Uj = ±a;i±a; 2 ±---±tu„-i, where k^ and u^ are the wave-vectors and the frequencies of the 
incident held. Phase-matching is used extensively to isolate specihc held-matter interaction 
pathways in non-linear spectroscopy. BB For example, in a degenerate 4-wave mixing 
(DFWM) experiment, [THTO] when the frequencies of all of the beams are the same, one 
commonly changes the directions of the incoming beams (noncollinear geometry) to spatially 
isolate the signal of interest. Although DFWM experiments provide a wealth of information 
on diherent resonant and non-resonant properties of matter, they cannot be used to study 
isolated microscopic systems whose dimensions are smaller than the wavelength of the light; 
the uncertainty in the momentum of the photons emanating from a microscopic object 
prohibits spatial separation of the different non-linear signals from each other and also from 
the excitation beams. This problem could be circumvented by using phase cycling schemes 
in optical spectroscopy, min] in analogy to multi-dimensional NMR spectroscopy. [12] In 
a typical phase cycling scheme, the relative phases between the different laser pulses are 
controlled by using pulse shapers and a set of measurements are done at different relative 
phases ( i.e., the relative phases in a single measurement are hxed). Addition of the spectra 
acquired at the well dehned set of relative phases selects the desired non-linear signal. aim 
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[H] The relative phases between the pulses can also be changed continuously (dynamic phase 
cycling), as it is done in the technique of phase modulation, [S], [IS] to separate the different 
signals arising from different held-matter interaction pathways. 

Phase cycling and phase modulation techniques have mainly been used in optical two- 
dimensional action spectroscopy. Ullsl na El HSl ETj In these experiments, one detects the 
action of four held-matter interactions on the incoherent signal, such as huorescence [Tl [5] fl^ 
Eiiisiini or photocurrent.[IHllin] However, these techniques can be used more generally in 
isolating different coherent and incoherent signals arising from the non-linear interaction of 
the laser pulses with matter. Motivated by the generality of the technique and its anticipated 
use in the spectroscopic investigation of isolated micro and nano-sized systems and devices, 
we present in detail the concepts and principles on the use of phase modulation in non-linear 
spectroscopy. Similar principles can also be applied in experiments that use phase cycling. 

The polarization P induced in a medium by an electromagnetic held E can be written 
as 

OO 

m = ( 1 ) 

i=l 

where Xi is the order susceptibility of the medium. Although a propagating electric 
held has both spatial and temporal dependence, we ignore the spatial dependence in the 
following derivation in order to simplify the notation. This simplihcation is justihed, as 
we use frequency rather than wavevector to separate the diherent non-linear signals in the 
phase modulation techniques. The total electric held due to a sequence of phase modulated 
pulses can be written as 

E{t) = ^ A(t)[exp{-i(a;j)t} -Fexp{i(a;j)t}], (2) 

j 

where loj = uj + cf)j, and (j)j is the modulation frequency of the pulse. In Eq.(|^, we have 
assumed that the amplitudes and the wavelengths of the diherent pulses (before the phase 
modulation) in the sequence are the same. [20] In a typical experiment, we use acousto-optic 
modulators (AOMs) to continuously sweep the relative phases of the optical helds. First- 
order dihraction from an AOM simply shifts the frequency of the incident laser beam, i.e. 
u u + cf). Here 0 is the frequency of the acoustic-wave in the AOM. Note that one can 
equivalently choose another hrst-order dihraction geometry with u — (p for all the beams 
to achieve the same results. Substitution of E{t) in Eq.([^ shows that the polarization 
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induced by n laser fields oscillates with frequencies icji ± 0 J 2 ± ^ oJn- It is possible to 

continuously sweep the optical phases by other means, such as reflection from a vibrating 
surface or an electro-optic modulator. However, such methods produce undesired sidebands 
at the multiples of the modulation frequency; uj ^ uj ± ncf), where n G {0, ±1, ±2, 

As we demonstrate below, this complicates the separation of different contributions to the 
non-linear signal. 

We provide examples that include sum-frequency generation, two-photon fluorescence, 
and two photon photocurrent to illustrate how phase modulation can be used to separate 
the various signal contributions to nonlinear spectroscopic experiments. 

II. PHASE MODULATED SUM FREQUENCY GENERATION 

In Figjg a), we show a schematic of the setup we used to investigate sum frequency 
generation by phase modulation. A pulsed laser beam is split into two beams, each of 
which passes through one of the arms of a Mach-Zehnder interferometer. Each arm of the 
interferometer has an AOM that modulates the phase of the beam. The beams are then 
recombined collinearly, and focused onto a non-linear crystal (NEC). The temporal overlap 
between the pulses is optimized by using a delay line (DL). The second harmonic and sum- 
frequency signals from the NEC are isolated from the fundamental beams by using a spectral 
filter (FIE), and detected by a photodiode (PD). The signal from PD is analyzed using a 
generalized lock-in amplifier fGEIAi. [2TlI23] 

The different second order non-linear processes induced in the NEC can be obtained from 
the expansion of the expression given in Eq.(|^: 

P 2 (t) oc {Ei(t) + E 2 (t)Y oc (exp{-i(ci; -F 0i)} exp{-i(ci; -F 02)} + c.c-Y 
(x2 + 2 [exp{-i (02 - 0 i)t} c.c] 2 [exp{-i( 2 a; + (pi + 02 )t} + c.c.] 

-f [exp{-i2(ci; 0i)t c.c.] -F [exp{-i2(ci; -F 02)t} + c.c.], (3) 

where c.c. stands for complex conjugate. The DC term and the terms with the frequency 
02 - 01 in Eq.([^ are due to optical rectification. Note that the optical rectification in 
conventional second harmonic generation (SHC) only has the DC component. When phase 
modulation is used, the optical rectification also produces an AC field at frequency 02 — 01- 
The terms with the frequency 2a;-|-0i 4-02 are due to sum frequency generation (SFC), while 
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FIG. 1. (a) Schematics of the setup used to analyze the sum frequency generation using phase 
modulation. OS: oscillator, OD FIL: variable optical density filter, BS: beam splitter, AOM: 
acousto-optic modulator, DL: delay line, M: mirror, NLC: non-linear crystal, FIL: filter, PD: 
photodiode, SIG: signal and GLIA: generalized lock-in amplifier. The AOMs in the two arms of 
the Mach-Zehnder interferometer modulate the phases of the beams by (j)i and (/> 2 , respectively. 
The two beams are recombined collinearly by BS2 and focused onto NLC. The second harmonic 
and sum frequency generated by NLC are spectrally filtered from the fundamental beams by FIL, 
and detected by PD. The signal from the PD is analyzed by GLIA. PD2 is used to record the 
intensity of the excitation beams, (b) Laser spectrum. 

the terms with the frequencies 2(a; -|- 0i) and 2(a; -f 02) are due to SHG from the individual 
beams. Only the contributions due to SFG and SHG in Eq.([^, which propagate as electric 
held Ep, are detected by the photodiode. The intensity of the signal, / oc recorded by 
the photodiode is given by Eq.Q: 

/ oc [exp{—z4(a; -|- 0i)t} + c.c.] -f [exp{—i4(a; -|- 02)0 + c.c.] 

3[exp{-i(4i:n 4- 20i -f 202)t} -F c.c.] 2[exp{-i(4a; 30i -F 02)0 + c.c.] 

-F 2[exp{-z(4a; + 0i + 302)0 + c-c.] + [exp{-z2(02 - 0i)t} + c.c.] 

-F 4[exp{-O02 - 0i)O + c.c.]. (4) 

Although Eq.(|^ has many oscillating terms, all of the terms that oscillate at high fre¬ 
quencies, i.e. the terms containing u, are hltered by the electronic circuit, so that only the 
terms with the frequencies 02i = 02 — 0i and 202i given in Eq.([^ are detected by the GLIA: 

Jdet oc [exp(-i202it) -1- C.C.] -F 4[exp(-i02if) + c.c.]. (5) 
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It is perhaps counterintuitive that the detected non-linear signals oscillate at the relative 
modulation frequency, 02i, and twice the relative modulation frequency, 2021- Although 
we have used a rather lengthy derivation to arrive at this relatively simple result, one can 
also use Feynman diagrams to reach the same conclusions. These diagrams also help us to 
rationalize why we observe the signals at the frequencies 02 i and 2021 - 

In Figs|^(i., ii.), we show the Feynman diagrams of SHG processes that lead to the oscilla¬ 
tion of the photocurrent from the PD at 2021- SHG processes described by pathways i. and ii. 
generate wave-packets \'ipa) ~ exp{—i2(a;-f0i)t} and \'ilJb) ~ exp{—f2(a;-|-02)O) respectively, 
in the PD. The interference between the two wave-packets ('0fc|'0a) generates photocurrent 
with phase signature exp(i2(;/)2it). Note that the conjugate wave-packet interference term 
('0a|'0b) produces a photo-current signal with phase signature exp(—^2021^). These two 
signals combine to produce a real-valued photocurrent at exp(i 202 it) -|- exp(—i 202 it) oc 
cos(202it). Other interference terms (t/’alV’a) and produce DG signals, which can be 

removed using a high-pass hlter preceding the GLIA. 
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FIG. 2. Feynman diagrams of SHG (i, ii) and SFG (iii) processes. The field due to i. creates a 
wave-packet and the field due to ii. creates a wave-packet Ik'b) in the detector. The interference 
between the two wave-packets, (k^blV’a)) leads to the oscillation of the signal at 2(^21- Similarly, 
field due to iii. creates a wave-packet \'ipc) in the detector. The interference of the wave-packets 
(V'alV’c) and (k’bl'fAc) generate the signal at (j)2i- 

In Figj^iii.), we show the Feynman diagram corresponding to the SFG process. Interfer¬ 
ence between the wave-packets generated by the SFG, I'll’c), and SHG, liJa), gives rise to the 
AG signal at 02i- Note that the interference of \'ipc) and I'lpb) also contributes to the signal 
at the frequency 02i- Thus, the signals at 2021 and 02i allow us to separate the second- 
order non-linear processes of SHG and SFG, respectively. We further note that optimal 
wave-packet interference due to the held-matter interactions of each beam requires that the 
beams are collinear. 
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We next demonstrate an experimental observation of the signals at frequencies (p 2 i and 
2021 that arise from SHG and SFG. We use a Synergy Ti:Sapphire oscillator from Femtolasers 
as the optical source in the setup [see FigB a)]. The oscillator delivers 10 fs pulses at a 
repetition rate of 80 MHz. As example laser spectrum is shown in Figj^b). The spectrum 
is centered ~780 nm and has a band-width of ~135 nm. AOMl is operated at 0i = 55 MHz 
and AOM2 is operated at 02 = 55.17 MHz. We used a 100 nm thick BBO crystal as the NLG. 
The SHG and SFG signals generated by the BBO crystal are hltered from the fundamental 
beams using a spectral hlter (86-973, FL hlter from Edmund optics). The hltered SHG 
and SFG signals at ~370 nm are then detected using a GaP photodiode (PDA25K-EG, 
Thorlabs). We note that linear absorption of the fundamental beams by a low band-gap 
photodiode would also produce a signal at the relative modulation frequency 02i. Therefore, 
a spectral hlter is necessary to separate the non-linear signal at the modulation frequency. 
In addition, as GaP does not absorb photons at the fundamental wavelength (700 - 950 
nm) the chance of detecting the linear signal becomes negligible. We also point out that 
the fundamental beams that leak through the hlter are too weak to induce any two photon 
absorption in GaP. 

In Figg a), we show the fast Fourier transform (FFT) of the SHG and SFG signals 
produced by the BBO crystal, and recorded by the photodiode. The FFT shows two peaks, 
one at 0.17 MHz, which is the relative modulation frequency of the AOMs (02i), and the 
other at 0.34 MHz (202i). The ratio of the amplitudes of the signals, S'i/S' 2 , is about 6:1, 
which is greater than 4:1, as predicted by Eq.(g. Two factors in the experimental conditions 
give rise to the deviation of the signal amplitudes from the ideal case given by Eq.(g. We 
have assumed that the amplitudes of the two beams, Ai and A 2 , are equal [Eq.(g] in 
deriving Eq.(g. When Ai ^ A 2 , the ratio S'i/S '2 > 4. This is also true when the two beams 
after the interferometer are not made to be perfectly collinear. In our experiments, Ai and 
A 2 are not exactly equal, as the beamsplitters used do not have 50/50 splitting ratio over 
the entire spectrum of the laser. It is also challenging to maintain the identical transverse 
beam prohles after the AOMs. Hence, perfect spatial overlap of the beams has not been 
achieved in our experiments. 

The dependence of the signal intensities at 0.17 MHz and 0.34 MHz on the intensity of the 
laser beams are shown in Figgb). Each data point was averaged over 20 measurements. The 
error bars in the data are the standard deviation of the different measurements. We did not 



FIG. 3. (a): FFT of the SHG and SFG signals recorded by the photodiode, (b): excitation intensity 
dependence of the signal at </> 2 i = 0.17 MHz and 2 (j)2i = 0.34 MHz. The points show the original 
data, and the lines are quadratic fits to the data points. 


directly measure the intensity of the laser beams in the measurements. Instead, we measured 
the amplitude of the modulated signal at 0.17 MHz using a silicon photodiode (labeled PD2 
in FigQ. This signal scales linearly with the intensity of the laser beams. A variable OD 
hlter (OD FIL in FigQ was used to change the intensity of the laser beams. The x-axis in 
Figj^b) is normalized by the maximum amplitude recorded by PD2, which corresponds to 
an average excitation power of about 10 mW (0.13 nJ per pulse) that was focused onto a 
spot size of ~ 20 /rm. The red points are the amplitudes of the non-linear signal recorded 
by PDl at 0.17 MHz, and the blue line is the quadratic equation f{x) = a + b{x — c)^ £t 
to the data points. Similarly, the data and the quadratic £t of the signal at 0.34 MHz are 
shown as black points and a green curve, respectively. The signal at 0.34 has been rescaled 
and shifted to fit in the plot. Both signals increase as a quadratic function of the laser beam 
intensity. These results show that the coherent second order non-linear signals induced by 
a pair of phase modulated collinear beams generate signals at 02 i and 202 i- 


III. PHASE MODULATED TWO-PHOTON FLUORESCENCE 

The strength of a two-photon absorptive transition by a fluorophore is given by 


S,p(t) X (Eiit) + E2(t))\ 


( 6 ) 
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Derivation of the full expression for S 2 P is similar to the derivation of P 2 - As in the case 
of SHG and SFG, we hnd that the intensity of the fluorescence induced by two-photon 
absorption is modulated at the frequencies 02i and 202i. In fact, the intensity of the de¬ 
tected two-photon fluorescence is also given by Eq.([^. Again, we use Feynman diagrams to 
rationalize this apparent similarity. 
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FIG. 4. Double sided Feynman diagrams of different pathways in two-photon fluorescence by phase 
modulation. Symmetric interaction of the flelds from the two beams (i) gives rise to the fluorescence 
at 2021. Asymmetric interactions of the flelds (three fields from one of the beams and one from 
the other) give rise to the fluorescence at 02 i. 


In Figj^ we show the Feynman diagrams of the different held interactions that give rise 
to two-photon fluorescence. The fluorescence intensity is modulated at the frequency 202i 
when one photon from each of the two beams is absorbed by the fluorophores, as depicted 
in FigQi). When three flelds from one of the beams, and the fourth held from the other 
beam participate in the absorption of two photons, the fluorescence is modulated at the 
frequency 02i. There are similarities between Fig|^i,ii) and FigQi). In the case depicted by 
Figgi ,ii), the two pathways generate two coherent second harmonic signals, which interfere 
in the detector to generate the photocurrent at 202i. For the process shown in FigQi), the 
two signals interfere in the sample itself to generate intensity modulated fluorescence at the 
frequency 202i. Similarly, the processes shown in Figj^i,iii) and FigQii), and Fig|^ii,iii) 
and FigQiii) produce analogous signals. 

In Fig^ a), we show the modulation frequencies of fluorescence signals resulting from a 
millimolar solution of Rhodamine 6G in water. The experimental setup is similar to that 
shown in FigQa). For this experiment, the NLG is replaced by a cuvette with the solution 
of Rhodamine 6G sample. The two photon fluorescence is collected at 90° to the optical 
axis of excitation, and focused onto an avalanche photodiode (APD) (LGSA3000-01, Laser 
Gomponents GmbH). A band-pass Alter at 550 nm (84-772, Edmund Optics) was placed 
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just before the APD to reject residual scattering of the fundamental beams. The FFT of the 
signal exhibits two sharp peaks at 0.17 and 0.34 MHZ, which correspond to the frequencies 
021 and 2021, respectively. The ratio of the amplitudes at the two frequencies is about 6:1. 
As for the case of the coherent signal from the BBO crystal, the ratio deviates from 4:1. 
Here too, the deviation can be attributed to the lack of perfect spatial overlap of the beams 
due to the differences in the beam prohles. 



Excitation Intensity [Arb. Units] 



FIG. 5. (a): FFT of the two-photon fluorescence from the solution of Rhodamine 6G, recorded 
by the photodiode, (b): Excitation intensity dependence of the signal at 02i =0.17 MHz and 
2021 = 0.34 MHz. The points show the original data, and the lines are quadratic fits to the data. 


In Figj^b), we show the amplitudes of the signals at 0.17 MHz (red points) and 0.34 
MHz (black points) as a function of excitation intensity. The blue and the red lines are the 
quadratic functions £t to the data. We see that the fluorescence at both the modulation 
frequencies increases quadratically with excitation intensity, which again demonstrates the 
nonlinear nature of these signals. Although these results are easily explained using Eq. 
they demonstrate an important application of phase modulation in non-linear spectroscopy. 
For example, only the signal at the frequency 202i has been used in multi-photon microscopy 
by phase modulation. m Our results and analysis show that the non-linear signal at 02i is 
more than four times stronger than the signal at 202i. Therefore, in microscopy where the 
non-linear signal can be spectrally isolated from the linear signal, one can also monitor the 
signal at the modulation frequency 02 i to achieve a better signal-to-noise ratio. 
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IV. PHASE MODULATED TWO-PHOTON PHOTOCURRENT 


We next present an example of the use of phase modulation in investigating non-linear 
optical effects in devices, such as photodiodes. As in the case of two-photon absorption 
induced fluorescence, when phase modulated beams interact with a two-photon photodiode, 
the resulting non-linear photocurrent is modulated at the frequencies <p 2 i and 2021 - In 
Figg a), we show the modulation frequencies of two-photon photocurrent from the GaP 
photodiode. These measurements were performed by replacing the NLC in the setup shown 
in Fig0 a) by the GaP photodiode. The signals at 0.17 and 0.34 MHz are similar to the 
signals obtained from SHG and two-photon fluorescence. We observed that the signals at 
both of these modulation frequencies increased quadratically with the excitation intensity 
[see Figgb)], which shows that both of the signals arose from the absorption of two photons 
by the photodiode. 



FIG. 6. (a): FFT of the two-photon photocurrent from GaP photodiode, (b): Excitation intensity 
dependence of the signal at 02i = 0.17 MHz and 2021 = 0.34 MHz. The points show the original 
data, and the lines are quadratic fits to the data. 


V. DISCUSSION AND CONCLUSIONS 

We have shown that modulation of the relative phase of two beams leads to the mod¬ 
ulation of the intensity of the coherent and incoherent non-linear signals at multiples of 
the phase modulation frequency, n02i. The amplitudes of the coherent signals due to the 
interaction of the system with three helds (second order non-linear process) are modulated 
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at the frequencies 02i and 2021- We have also shown that the amplitudes of the incoherent 
signals due to the interaction of the system with four helds are also modulated at 02 i and 
2021- Theoretical analysis shows that the amplitude of the signal with modulation frequency 
021 is four times larger than the amplitude of the signal with the frquency 202i. In our ex¬ 
periments, we observe that the signal at 02 i is roughly six times larger than the signal at 
2021- We attribute this discrepancy to the imperfect overlap of the transverse beam prohles 
at the detector ( for the case of the coherent signals) and at the sample (for the case of the 
incoherent signals). Imperfect overlap occurs when the beams are not exactly collinear, and 
also when the modes of the beams do not exactly match with each other. 

Extension of the analysis of the held-matter interactions to higher orders reveals that the 
amplitudes of the coherent signals, which are due to n — 1 held-matter interactions, and 
the incoherent signals, which are due to n held-matter interactions, are modulated at the 
frequencies given by x<p 2 i, where x = n/2, n/2 — 1,..., 1, and n is an even integer. In the three 
examples presented in this work, we have not observed intensity modulation at frequencies 
higher than 202i, and the signals scale quadratically with the excitation intensity. We 
have thus been able to observe only the interactions that involve at most four held-matter 
interactions. This limitation is due to the maximum intensity of the excitation beams 
available from our source laser oscillator. However, laser beams from amplihed systems 
should be able to drive even higher order non-linear processes. The phase modulation 
technique could be useful to the investigation of diherent interaction pathways in such high 
order processes. Moreover, as the technique can be used to investigate the diherent non¬ 
linear contributions to incoherent signals, the method could also be used to investigate the 
spectroscopic properties of isolated systems, such as individual molecules, nano-particles, 
and devices. 
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